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Abstract 

A theory for the spin wave eigenmodes of a Dzyaloshinskii domain wall is presented. These walls 
are Neel-type domain walls that can appear in perpendicularly-magnetized ultrathin ferromagnets 
in the presence of a sizeable Dzyaloshinskii-Moriya interaction. The mode frequencies for spin 
waves propagating parallel and perpendicular to the domain wall are computed using a continuum 
approximation. In contrast to Bloch-type walls, it is found that the spin wave potential associated 
with Dzyaloshinskii domain walls is not reflectionless, which leads to a finite scattering cross-section 
for interactions between spin waves and domain walls. A gap produced by the Dzyaloshinskii 
interaction emerges, and consequences for spin wave driven domain wall motion and band structures 
arising from periodic wall arrays are discussed. 
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I. INTRODUCTION 


The Dzyaloshinskii-Moriya interaction (DMI) is an antisymmetric contribution to the 
exchange energy that can exist in spin systems that lack inversion symmetry. [IH3] Spin 
textures stabilized by DMI are of special interest due to the possibility of new technological 
applications. @HS] For room temperature operation, interface DMI is of particular importance 
in terms of thin film structures based on transition metals, compatible with traditional 
spintronic devices. It has been shown in perpendicular materials that because of DMI 
compensation of the dipole-dipole interaction at the center of a domain wall, a Neel-type 
domain wall is favored with important enhancements of domain wall stability and mobility. [7, 

i 

Small amplitude fluctuations of the magnetisation about equilibrium- spin waves- in 
systems with DMI have been studied experimentally, [TO] and theoretically [lTJ [12] for 
interface DMI films. A key feature is nonreciprocity of frequency as a function of propagation 
direction for finite wavelength spin waves, which emerges from the chiral symmetry breaking 
DMI. In the present work we examine the spin wave dispersion in a film containing a DMI 
stabilized Neel wall. We find that spin waves are partially reflected by the domain wall 
due to an extra chiral term in the potential associated with the domain wall. Local modes 
are found, and we show that DMI drives a hybridisation of traveling spin waves with these 
localized states. 

The hybridisation produces an energy-split dispersion with a gap magnitude that is pro¬ 
portional to the magnitude of the DMI. We illustrate this effect with a suggestion for a 
magnonic crystal produced by a periodic array of domain walls in a nanowire with DMI. 
Without DMI the walls are Bloch-type walls which are known to represent reflectionless 
potentials for spin waves traveling through them. [IMS] This results in a gapless band 
structure. However, with DMI the stable configuration becomes an array of Neel-type walls 
with a modified potential for the spin waves. The modified potential is no longer reflection¬ 
less and standing waves appear at the edges of the first Brillouin zone producing gaps in the 
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FIG. 1. (Color online) Geometry considered for the Neel-type Dzyaloshinskii domain wall. Xq 
denotes the position of the wall center along the x axis. Translational invariance is assumed along 
the y direction and the magnetization is taken to be uniform across the thickness of the film in the 
z direction. 

band structure. 

This article is organized as follows. In Section II, the model and calculations involving 
the static domain wall profile are presented. In Section III, the spin wave eigenmodes of 
the Dzyaloshinskii domain wall (DDW) are computed using a variational method in the 
continuum approximation. Consequences of these eigenmodes are then explored in Section 
IV, where the reflection and transmission coefficients for propagating spin waves through 
the domain wall are computed and band gaps in associated with periodic wall arrays are 
discussed. Finally, a discussion and some concluding remarks are given in Section V. 


II. MODEL AND STATIC WALL PROFILE 

An ultrathin ferromagnetic wire is considered in which a domain wall separates two 
uniformly-magnetized domains along the x axis, as shown in Figure ([Tj) . The uniaxial 
anisotropy, K u , is taken to lie along the z axis, perpendicular to the film plane, while a 
transverse anistropy resulting from volume dipolar charges, K±, is present along the x axis. 
In addition, an interfacial Dzyaloshinskii-Moriya interaction D is also considered, with a 
form consistent with a multilayered system with a heavy-metal subtrate. mm The form 
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of this interaction can be written in terms of the Lifshitz invariants L% = — rrij^p^ 

as D (L* x + L y zy ). The magnetization orientation, represented by the unit vector m, is 
parametrized using spherical coordinates as m = (sin 6 cos 0, sin 6 sin 0, cos 9), where 6 = 
9(r,t) and 0 = 0(r,t). The total magnetic energy of this system is given by the functional 
U[9( r),0(r)], 


U = dV [A((W) 2 + sin 2 0(V0) 2 ) 


+ (K u + K l cos 2 0) sin 2 6 

fd9 , do . , 

+ V [ — cos 0 + — sill 0 
\ox ay 


d(j) 


dcj) . 


+ - sin 26 —— cos 0 — -- sin 


dy 


dx 


, ( 1 ) 


where A is the exchange constant. The static profile of the domain wall is determined by 
the solution to the Euler-Lagrange equations associated with the functional in Equation ([Tj) , 
which are obtained by setting the first-order functional derivatives to zero. By neglecting 
variations in the y direction and assuming the solution 0(x) = 0o, the nonlinear differential 
equations satisfied by the static wall profile ( 9 0 , 0 O ) are given by 


S—' S ‘ 

5U 


- - (K u + K ± cos 2 (0 o )) sin(26> 0 ) = 0, 

.dOr 


50 


—- = 0 =>• sin(0 o ) siir(6 , 0 ) ( K± cos(0 o ) + h— = 0. 


dx 


( 2 ) 

( 3 ) 


Note that the second equation above is satisfied by the 0o ansatz for finite values of the 
DMI, D ^ 0, only if the domain wall assumes a pure Neel profile (0 O = 0,7r). By assuming 
a Neel wall state, the solution for 0 o (x) to be written as 

x - X 0 


0q(x) = 2 tan 


-i 


exp 


A 


( 4 ) 


where A = \JAf (K u + K ±) is the domain wall width parameter and Xq denotes the wall 
center. The solution with the positive sign in the argument of the exponential function gives 
the configuration illustrated in Figure ([Tj) . With this solution, the total domain wall energy 
(JT|) can be evaluated to be 


= U[6 0 , 0 O ] = Ay/A (K u + K ± ) =F t tD, 


( 5 ) 
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where the negative sign corresponds to the solution (j) 0 = n and the positive sign to 0o = 
0, which indicates that left-handed Neel walls are preferred energetically for D > D c = 
AXK±/n > 0. 


III. SPIN WAVE HAMILTONIAN 

The magnetic energy functional can be expanded up to second order in small fluctuations 
(56, 5<f>) around the stable configuration (6 0 , cp 0 ) to obtain the spin wave Hamiltonian. [ 131 [17] 


K 


5H = — / dx 56 Vp(x) 56 + 5(f> \Vp(x) — A(x) — n] 5<p. 


K 


( 6 ) 


The energy of the fluctuations is described by the operators Vp(x) = [—A 2 <9 2 + 1 — 
2sech 2 (x/A)], A(x) = (D k/XK±) sech(x/A) and k = K±/(K U + K±). The Schrodinger 
type operator, Vp(x), has been widely studied and is used to describe spin waves in a Bloch 
type domain wall. [131 H3J Solutions to these operator include a single bound state, 

€ loc ( x ) = sec h(x/A), (7) 

with zero corresponding energy, and continuum-traveling states, 

£k(x) = : e lkx [tanh ( x/X ) — ikX], 


( 8 ) 


with eigenenergy given by ujk = 1 + k 2 A 2 . The above states form a complete orthonormal 
set, 

[ dV CkCioc = 0 , 

J r (9) 

y dV e k U = 5 Km . 

From Equation (J6]) it can be deduced that k introduces a constant ellipticity in the precession 
of the fluctuations, and DMI introduces a spatially dependent one through the A (A) term 
so that it is not trivial to find a basis that diagonalizes the spin wave Hamiltonian. We 
propose a linear superposition of the local and traveling modes, 


x(x) = 5(p(x) + i56(x) = ici 0 S 0C (x) + ^ d k ^ k (x), 


( 10 ) 
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to calculate the spin wave energy. After the space integrals are computed, we find 

SH = -cL (^ + A' ± ) + £ [A' k dld k + B' k (d‘ k d‘_ k + d t d_ k 

+CfcQoc(dfc + d* k )} + Ukm d* k d m + V km (d k d m + d k d^), (11) 


where the coefficients are given by 


km 


I<i 


A' k = u k (K u + K ± ) - 

Cfc = I dV ZiocA(x) 4, 
Ukm = / dV £lA(x)£m, 


( 12 ) 


^fcm= / dV^ k A(x)U+C k ^(x)0- 

The A' fc and terms denote elliptical spin precession as a result of the transverse anisotropy 
and correspond to the usual terms found in the Bloch wall case.pZQ E] The C k , U km and 
V km terms are proportional to the strength of D and depend on k because these terms result 
from the spatial dependent cllipticity. C k represents the coupling between the local and the 
traveling modes, it is small compared to the other terms so it will not be considered. U km 
and V km are scattering terms that describe the transition from a state with momentum hk 
to another state with hm. If we focus on the maximum scattering strength then the specific 
form of the coefficients, U km ~ sech (k — m) and V km ~ sech(/c+m), allows us to approximate 
U km and V km by delta functions 5 km , S k - m respectively. We can then approximate 5H as 

5H = -c 2 oc + Kj^j + ^ A k d* k d k + B k (d* k d*_ k + d k d_ k ). (13) 


The first term on the right hand side of Equation (13) can be related to the domain wall 


mass by p 2 /(2rri]y ) ~ cf oc (qy + K±)i where = l/[2(^ + Kj_)\ is the Neel-type domain 

wall mass. [20, [2T|. The mass in a Bloch-type wall is mg = l/(2A'j_) so < tub which 

agrees with a higher mobility in DDWs.[6j The rest of the coefficients are 

K ± 7 tD (1 + 2k 2 X 2 ) 

(14) 


A k cu k (A u 
K, 


K i ) - 


Bu — 


2 4A 
7 tD (l + 2k 2 X 2 ) 


u k 


8A 


u k 
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This Hamiltonian can be diagonalized by means of a Bogoliubov transformation, c k = u k d k — 


u k dh u k = V (-^fc ± to obtain 

6H = -<Z'(^ + K J )+'£,m<iici', (15) 

' ' k 

where the frequency Q k is given by 

„ (K u + K ± )a 3 I 7 ^ 75 (1 + 2i 2 A 2 

*“ ft K 2\(K u + K ± ) (1 + ftW) 

where a ~ 0.3 nm is the lattice constant, ft is now possible to explicitly write the spin waves 
eigenmodes in terms of the amplitudes Ck, c*_ k and the local and traveling modes 



x{x) = icioS oc (x ) + c k ul + c*_ k u k )£ k (x) 

k 


(17) 


IV. BAND STRUCTURE IN PERIODIC WALL ARRAYS 


The scattering potential for spin waves in a Bloch (D = 0) domain wall is represented by 
Vp(x) which is reflectionless but leads to a a phase shift when spin waves propagate through 
it- H51 Vp(x) corresponds to a specihc case of the so called modified Poschl-Teller 
Hamiltonian, [23] 

[— ot 2 d1 — 1(1 — 1) sech 2 (ic/a)] if) = eip. (18) 


The parameter l describes the depth of the potential well, a has units of distance and e 
is a dimensionless energy. For a Bloch-type wall, l = 2 and a = A. The transmission 
and reflection coefficients related to the wave propagation across this potential have been 
calculated for this Hamiltonian as a function of the depth 


l«p = 


I T? = 


P 


(19) 


1 + P 2 ’ 1 + P 2 ’ 

with p = sinh(7rfco;)/sin(7rZ).[24] From this result it can be seen by inspection that for l e N, 
\R\ 2 is zero. For the Dzyaloshinskii domain walls, the Hamiltonian is 

D sech(x/A) 


A 2 cr — 2seclr(ic/A) 


A (K u + K ± )_ 


X(x) = Ex(x), 


( 20 ) 
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where the dimensionless energy is E — + k — 1. The D term modifies the depth 

of the potential but not its form, ft possible then to relate the parameter l with D , 


i = i 
2 


1 + W 1 + 4 ( 2 + 


D 


X(K U + K ± ) 


( 21 ) 


Two transmission coefficients were calculated as a function of the wave vector k using Equa¬ 
tion (21) for different values of D and are shown in Figure ([2]) along with numerical sim¬ 
ulations to verify our theory. The numerical calculations were performed within a micro- 
magnetic model. The calculations were done with the code mumax3.[25lj The standard code 
includes the interface DMI term but was modified to include at the same time the in-plane 
and out-of-plane anisotropies. The parameters used were A = 16 pj/m K\ = 18 kJ/m 3 , 
K u = 0.5 MJ/m 3 and A = 5.55 nm. The system was discretized in cells of 1.5625 x 1.5625 x 1 
uni 3 . The geometry coincides with the one showed in Figure 1 and the system size was 
12800 x 50 x 1 nm 3 with periodic boundary condition in y direction. To simplify the 
analysis and comparison with the analytical model the calculations were performed without 
damping term and demagnetizing field. A domain wall was introduced at the center of the 
sample and then the system was excited with a monochromatic point source of 50 mT ap¬ 
plied field, 1950 nm away from the domain wall. The amplitudes were calculated comparing 
the average envelope of the spin waves at both sides of the domain wall at the initial stages 
of the propagation. As D increases significant reflection is found for larger values of k. This 


is a direct result of the scattering terms in Equation (11). 


As a result of the DMI, the scattering potential associated with the domain wall produces 
reflection in the spin waves propagating through it. As such, the band structure for a lattice 
of DDWs presents gaps at the edges of the Brillouin zone because of Bragg reflection, which 
is not present for Bloch-type walls for which no reflection occurs. To see this we consider 


a periodic array of DDWs and Fourier transform Equation (20) using also Bloch’s theorem 
on x( x ) to obtain the central Equation 


({K u + K x )X 2 k 2 - E ) C(k) + U G C(k -G)= 0, 

G 
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( 22 ) 











FIG. 2. (Color online) Transmission coefficient for D = 1.3 mJ/m 2 (red) and D = 2.6 mJ/m 2 


(black). The solid lines result from using equation (21) and the points are numerical simulations 


where Uq are the Fourier coefficients of the potential. (26] The period of the DDW crystal 
can be determined with the Kooy-Enz formula that describes the stray field energy for an 
arrangement of parallel band domains separated by domain walls of zero width. [2T] For a 
particular case of D — 2.6 mJ/m 2 a minimum him thickness of approximately 2 nm is found 
with a period of L — 100 nm. There is a compromise between the him thickness and the 
period, since the minimum him thickness and period increase as D decreases. 


Equation (22) represents an infinite set of equations connecting the coefficients C{k — G) 
for all reciprocal lattice vectors G. These equations are consistent if the determinant of the 
coefficients is zero. It is often only necessary to consider the determinant of a few coefficients. 
For our calculations an 11 x 11 matrix is used to numerically solve the central equation. 
The calculated band structure of domain wall crystals is shown in hgure [3j Gaps in the 
band structure are a consequence of Bragg reflection and a direct result of the DMI. Figure 
0 shows the hrst gap frequencies as a function of D in k = 0 and k = tt/L. 


V. DISCUSSION AND CONCLUDING REMARKS 


Rehection of spin waves by a domain wall is found when the interface form of DMI is 
included. It is a result of the stabilization of a Neel wall as the stable conhguration, and 
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FIG. 3. (Color online)Band structures of a domain wall crystal, (a) Bloch wall (D = 0).(b) Neel 
wall (D = 1.56 mJ/m 2 ). 



FIG. 4. (Color online) Frequency gaps A F at the Brillouin zone boundary as a function of D c . 
L = 100 nm' 1 is the period of the crystal. 

of an extra chiral term in the Hamiltonian that changes the Poschl-Tcller potential and 
scatters the spin waves. Reflection results in energy gaps in the band structure of a periodic 
array of domain similar to the ones found in a magnonic crystal. Our proposed model offers 
an alternative method for constructing a magnonic crystal without the need to build the 
metamaterial although we recognize the difficulty of stabilizing the domains. Moreover, the 
gaps only depend on D so there is only one parameter to control. It is worth noting that 
the bulk form of DMI favors a Bloch-type wall configuration and no extra term is found in 
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the spin wave energy. This last statement agrees with previous claims that the reflectionless 
feature is very robust. [28] 
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